Chemical-potential route for multicomponent fluids 
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The chemical potentials of multicomponent fluids are derived in terms of the pair correlation 
functions for arbitrary number of components, interaction potentials, and dimensionality. The 
formally exact result is particularized to hard-sphere mixtures with zero or positive nonadditivity. 
As a simple application, the chemical potentials of three-dimensional additive hard-sphere mixtures 
are derived from the Percus-Yevick theory and the associated equation of state is obtained. This 
Percus-Yevick chemical-route equation of state is shown to be more accurate than the virial equation 
of state. An interpolation between the chemical-potential and compressibility routes exhibits a better 
performance than the well-known Boubh'k-Mansoori-Carnahan-Starling-Leland equation of state. 
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I. INTRODUCTION 

It is well known that in a fluid made of particles 
interacting via a pair-wise potential, the most rele- 
vant statistical-mechanical quantity in equilibrium is the 
pair correlation function, or radial distribution function 
(RDF), g(r) 0, It is defined as the average number 
density at a distance r from a certain reference parti- 
cle, relative to the global density. Apart from account- 
ing for the structural correlation properties of the fluid, 
the RDF g(r) allows one to obtain a number of thermo- 
dynamic quantities {ipi} by means of elegant formulas. 
Since those thermodynamic quantities are connected by 
differential relations, one would expect to get the same 
macroscopic description, i.e., the same free energy A, re- 
gardless of the specific route g(r) — > ipi — > A followed 
[3j. However, this is not necessarily the case when an ap- 
proximate RDF is used as a starting point, what results 
in the well-known thermodynamic inconsistency problem 

a- 

To be more specific, let us assume an s-component sys- 
tem made of N = N a particles (N a being the num- 
ber of particles of species a) enclosed in a d-dimensional 
volume V. The total number density is p = N/V and 
the partial number densities are p a = N a /V = x a p, 
x a = N a /N being the mole fractions. The pair inter- 
action potential, not necessarily isotropic, between two 
particles of species a and 7 is </> Q7 (r) = 7Q (— r). In 
terms of the pair correlation function of species a and 7, 
g ay (r), it is possible to express the pressure p as Q 
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Ca^y dr3 a7 (r)r-V0 Q7 (r), (1.1) 



where Z is the compressibility factor and (3 = 1/ksT (fcs 
and T being the Boltzmann constant and the tempera- 
ture, respectively). Analogously, the internal energy per 
particle u = U/N can be expressed as 



(r). (1.2) 



Equations (|1.1[) and (|1.2[) are the virial (or pressure) and 
energy routes to thermodynamics, respectively. The free 
energy A = Na (where a is the free energy per particle) 
can be (partially) obtained Q either from Eq. (|1.1[) or 
from Eq. (( 1 . 2|) by inversion of the thermodynamic rela- 
tions 



Z = p 
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Apart from the virial and energy routes, a popular 
route is the compressibility one. It reads 



a, 7 



-1 

Q7 



(1.5) 



where the element h ai of the matrix h is proportional to 
the zero wavenumber limit of the Fourier transform of the 
total correlation function /i Q7 (r) = g ai (v) — 1, namely 



(1.6) 



— P\J x a Xy I dr h a ~y (r) . 
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Combination of Eqs. (|1.3[) and (|1.5j) allows one to get the 
free energy from the RDF via the compressibility route. 
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Equations (jl.3[) and (JTT4J) are related to the derivatives 
of the free energy per particle a with respect to the total 
number density and to the temperature, respectively, but 
they ignore the composition-dependence of a. In fact, the 
partial derivatives of a with respect to the partial number 
densities are related to the chemical potentials: 



dpa s 



(1.7) 



While in Eq. (|1.3|) the free energy per particle a is seen 
as a function of the total number density p and the s — 
1 independent mole fractions {x±, X2, ■ ■ ■ ,x s -±}, in Eq. 
(|1.7[) a is seen as a function of the s partial densities 
{pi, p2, ■ ■ ■ , p s }- Adopting the former point of view, Eq. 
(|1.7[) can be rewritten as 

(1.8) 

This allows us to get the useful identity 



E 



( dpa \ 

V d P )fl,{Xa} 



(1.9) 



Recognizing that [9(pa)/9p]^ r , = a+p/p is the Gibbs 
free energy per particle, Eq. (|1.9[) is not but the funda- 
mental equation of thermodynamics 

Cross partial derivatives of the free energy should be 
independent of the order of derivation, which yields the 
following Maxwell relations involving the chemical poten- 
tial: 



V 9(3 ) {pa} ~ \d P J 0!{ } 



(1.10) 
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(1.12) 



(1.13) 



In Eqs. (fTTTT) and (fl~T2|) the replacement p a -> A Q /V 
must be done in order to express the chemical potential 
as a function of {N a } and V. 

The interesting question is, how can one obtain the 
chemical potential p v (and hence the free energy) from 
the knowledge of the RDF? This question has been ad- 
dressed in several textbooks [H, Q and classical papers 
@, 0, but the most general formula (valid for any di- 
mensionality, interaction potential, and coupling proto- 
col) seems not to have been derived yet. The first aim 



of this paper is the derivation of that chemical-potential 
route to thermodynamics [see Eq. (|2.23[) below] . The sec- 
ond goal is the application of the route to a mixture of 
(additive) hard spheres (HS) in the context of the Percus- 
Yevick (PY) approximation [see Eq. (|5.4[) below]. This 
extends to multicomponent systems, the work recently 
reported in Ref. @ . 

This paper is organized as follows. The formal deriva- 
tion of the chemical-potential route is carried out in Sec. 
HU As a simple test, it is checked in Sec. [Hi] that the 
use of the exact pair correlation function to first order 
in density provides the exact second and third virial co- 
efficients. Section IIVI is devoted to the specialization of 
the chemical-potential route to d-dimensional HS mix- 
tures with zero or positive nonadditivity. In the three- 
dimensional case, the knowledge of the contact values of 
the radial distribution function for additive HS mixtures 
in the scaled-particle theory (SPT) and PY approxima- 
tions is exploited to obtain the chemical potential. It is 
observed that, while the SPT yields a result thermody- 
namically consistent with the virial route, this is not the 
case of the PY result. Finally, the paper ends with some 
concluding remarks in Sec. I VII 



II. THE CHEMICAL-POTENTIAL ROUTE 

As stated before, we consider an s-component mixture 
with N a (a — 1, . . . , s) particles of species a and N = 
Y^a—i Aa total number of particles in a volume V. We 
will employ the short-hand notations N = {N\, . . . ,N S }, 
r N = {ri, . . . , rjy}, and dr^ = dri • • • dr^, r.j being the 
spatial coordinates of particle i. If we denote by $N(r Ar ) 
the total potential energy, the (canonical-ensemble) con- 
figurational probability density is @, Q 



QnGS.V) 



(2.1) 



where 



(P,V) = V- N /dr Jv e-^ r ") (2.2) 



is the configurational integral. The average number of 
pairs of particles (per unit volume) of species a and 7 
located at r a and r 7 , respectively, is 



^07 



> Q ,r 7 )= / dr Ar / 9 N (r 7V )^(S eiiCe (5 ej 



xS{ri - r Q )<5(r J - r 7 ), 



(2.3) 



where ej denotes the species particle i belongs to. We 
define the pair correlation function as g a7 (r a , r 7 ) = 
n Q7 (r Q , v 1 )jp a p 1 . Inserting Eq. (|2.ip into Eq. (|2.3p . one 
obtains 



T/-(JV-2) 

?Q7 (r Q ,r 7 ) = ^— ldr N e-^^ 



X(5(ri - r Q )S(r 2 - r 7 ), 



(2.4) 
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where, without loss of generality, particles i = 1 and j = 
2 are assumed to belong to species a and 7, respectively. 

Contact with thermodynamics is made through the 
free energy An of the system: 

A N (p,V) = -k B T\nZ N {p,V), (2.5) 

where the partition function i? N factorizes as 

Z n (/3,V)=Z$(i3,V)Qk(/3,V), (2.6) 



^zid 



iP,V) 



V 



N 



dN„ 



(2.7) 



being the ideal-gas partition function. In Eq. (|2.7|l , A Q = 
/i/V27rm Q fcsT (where /i is the Planck constant and m a 
is the mass of a particle of species a) is the thermal de 
Broglie wavelength. 

Let us now focus on a given species v. The associated 
chemical potential is 



_ dA N (p,V) _ id 

llv - dN v - ^ 



where 



id 



dN u [P " ' 



(2. 



(2.9) 



^7 = - 



= In 



dhiQ N (/3, V) 
dN v 
Qn(P,V) 
Qn+i(W 



(2.10) 



In the second line of Eq. (|2.10[) we have taken into ac- 
count that N v 1 and it is understood that the sub- 
script N + 1 means {Ni , . . . , N v + 1, . . . , N s }. Obviously, 
Qn+i((3 7 V) is given by Eq. (|2.2p with the replacements 
N -> N + 1 and N — > N + 1. Moreover, without loss 
of generality, we will assign the label i = to the extra 
particle of species u, so that dr^" 1 " 1 = drodr^. 

Now we assume that the potential energy is pair-wise 
additive, namely 



*n(0 = 



N-l N 

E E 



(2.11) 



$N+i(r 



N+l\ 



N 

E 

3=1 



^.(i-ci^ + SnK), (2.12) 



where ^ a7 (r a ,r 7 ) = Q7 (r 7 — r Q ) is the interaction po- 
tential of two particles of species a and 7. 

In order to establish a relationship between the chemi- 
cal potential and the pair correlation functions, it is con- 
venient to introduce a coupling parameter £ such that its 
value < £ < 1 controls the strength of the interaction of 



particle i — to the rest of particles. A similar charging 
process was employed by Onsager in a different context 
[H|. It is also closely related to the so-called Widom 
insertion method [ll| . 

In our specific problem, the interaction potential be- 
tween particles i = and j > 1 is ($1 (tq, Yj), with the 
boundary conditions 



0, e = o, 

<?W,(ro,r,-), f=l. 



(2.13) 



The associated total potential energy and configuration 
integral are 



N 



^(r^ 1 ) = (ro.r,-) + $ N (r w ), (2.14) 



N+1 (/?,F) = / dr^e"^^^. 

(2.15) 

Note that Eqs. ([2~14|) and (|2~T5]) (with < £ < 1) define 
a system of A^ + 1 particles and s + 1 species, one of the 
species being made by the single particle i = 0. Only in 
the limits £ = and £ = 1 does one recover a number of 
species s. Analogously to Eq. (|2.4|) . the pair correlation 
function of particle i = and any particle of species a is 

Si/a ( r i / ) r a) 



X(5(r - rv)(5(ri - r a ). 



dr W+l e -/i*^ 1 (r" + 1 ) 



(2.16) 



Obviously, 



$ «) f^J^). ? = ° f217) 

$N + i(r } - u N+1 (^ +1 ), e = i, (2 - 17) 



(2.18) 



ff«W«) = I 1 ' , v ! ? (2-19) 
[gwiiV,]^), £ = 1. 

Assuming that the potentials are differentiable 
with respect to £, the second line of Eq. (|2.10[) can be 
rewritten as 

/3/C = -/ d^lnQ^CW (2.20) 



Now, from Eqs. (|2.14|l and (|2.15[) we obtain 

N 

= -^-^^ / , Q / dr JV + 1 e-^*»^^ +1 ) 
xS^CrcrO, (2.21) 
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where again particle j = 1 is assumed (without loss of 
generality) to belong to species a. Making use of Eq. 

(me}, 



where f${r„,T a ) = e'^^*"'*^ - 1 is the Mayer func- 
tion @,[|. 

To first order in density. 



0«l»Q&i = -£l>" / dr " /dr Q ^(r„r Q ) y^(r v ,v a ) = 1 + P, j dr 7 /«> (r„, r 7 )/ 7Q (r 7 , r Q ) 



xa e 4«j(r„r Q ). (2.22) 
Combination of Eqs. (|2.20[) and (|2.22j) yields the desired 



result 



^ = Eft.f« / dr^W^W, (2.23) 

where we have taken into account the translational in- 
variance property <^,a(r„,r a ) = </>ia(r Q - r„). In terms 
of the cavity function y$(r) = ^(rje^to, Eq. (12331 
becomes 



d£ / dry^r) 



(2.24) 



Equation (|2.23p or, equivalently, Eq. (|2.24[) consti- 
tutes the chemical-potential route to thermodynamics. 
It expresses the chemical potential of species v in an s- 
componcnt mixture (the "solvent") in terms of the pair 
correlations of an (s + l)-componcnt mixture where the 
extra component consists of a single particle (the "so- 
lute") whose interaction with the rest of the particles 
is controlled by a coupling parameter £, which switches 
from £ = (the solute does not feel the presence of the 
solvent particles) to £ = 1 (the solute is indistinguish- 
able from a solvent particle of species v). The result is 
independent of the number of species s, the dimensional- 
ity of the system d, the pair interaction functions <j> ay (r) 

(which can be isotropic or not), and the protocol <01(v) 
followed to go from </4a(r) = to (01 (r) — </>„ Q (r). If 
the exact gi^2{r) is employed, one gets the exact chem- 
ical potential p v (and, hence, the exact free energy A) 
regardless of the protocol. This is illustrated in Sec. IIIII 
at the level of the third virial coefficient. On the other 
hand, if an approximate function g{y& (r) is used instead, 
not only the resulting approximate free energy will be dif- 
ferent from the one derived from any of the other routes 
(e.g., virial, energy, and compressibility), but it will, in 
general, depend on the choice of the protocol. Thus, the 
chemical-potential route adds an extra source of thermo- 
dynamic inconsistency. 



III. LOW-DENSITY REGIME 



The aim of this section is to exploit Eq. (|2.24[) when 
the pair correlation function is truncated to first order in 
density. First, let us rewrite Eq. (|2.24[) as 

ptf = -V- 1 j d^PcJ dr„ J dr a y$(r„,r„) 

x0 e /tt>(r„,r a ), (3.1) 



+o( P 2 ), 



(3-2) 



where / 7a (r 7 ,r Q ) = e 04>7<»( r 7> r <*) — l. Insertion into Eq. 
(|3~Tj) yields 

Pp°* = -V^ 1 ^2 p a J dr„ j &r a f va (r u ,r a ) 

a 

-V- 1 / d£/«)+0(A (3-3) 
Jo 

where we have called 

IP =J2p*P~f / dr„ / dr Q / dr 7 /$(r 

a, 7 

xd ( f^(v V) r a ). (3.4) 

Exchanging a <-> 7 and r Q O r 7 in the summation and 
the integral, ijp can be rewritten as 

IP = r pap i j AVu [ dVa /dr 7 / 7Q (r 7 ,r Q ) 

a, 7 J J J 

xdt[f<g(T Vi T a )f<$(r v ,r y j\ . (3.5) 
Therefore, Eq. ()3.3[) becomes 



a, 7 



where 



T va = / drf ua (r), 



(3.7) 



Q vai = J dr / dr'/ 7a (r - r')f va (r)f vl (r'), (3.8) 

and the translational invariance property has been ap- 
plied again. Note that both T va and Q vart are invariant 
under any permutation of indices. This guarantees that 
Eq. (fl~T3"l) is verified. 

Equations (|3.6|) - (|3.8|) show that, as expected, the spe- 
cific protocol does not play any role in the final result. 
From Eq. (|1.7p we finally get 



(3pa CX = -^^PaPyFay PaP 7 PuGa 7 u + 0(p 4 ). 

(3.9) 



ct,7 ct,7,i/ 



Making use of Eq. (|1.3I) it is straightforward to get 

Z = l + B 2 p + B 3 p 2 + 0{p 3 ), (3.10) 
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where the second and third virial coefficients are 

B 2 = -^-'^x a x 1 F ai , (3-11) 



where 



P^' 1 = d2 d v d da 0a a$ a l yoa((To a ), (4.6) 



(3.12) PiC' n = d2 d v d J2 Pa- / ddo^oVf/Oat^). (4.7) 



a, 7,1/ 



As expected, these are the exact expressions Q. 



IV. HARD-SPHERE MIXTURES 

Let us now particularize the chemical-potential route 
(|2.24[) to HS mixtures. In that case, 



" (r) = 6(r - Very), 



(4-1) 



where r = |r|, Q(x) is Heaviside's step function, and a aj 
is the range of the infinitely repulsive interaction between 
particles of species a and 7. 

We now need to introduce the extra particle i = 
(the solute) coupled to the remaining N particles through 
a coupling parameter < £ < 1 via the set of in- 
teraction potentials (r). According to Eq. (|2.13[) . 
cxp[-/3^(r)] = 1 and exp[-^£ri (r)] = 0(r - ov Q ), 
but otherwise a certain freedom to fix the protocol <j>^2 (r) 
exists. The most natural choice is a HS form, i.e., 



e-^W=e(r- ff tg), 
with (ria' = and ctIq = <r VOL . Therefore, 



9£ ' 



so that Eq. (|2.24j) becomes 

/3^ x = d2 d v d y^jQq / dcr 0ce o-q'^qq^oq), 



(4.2) 



(4-3) 



(4-4) 



where we have taken into account that the integral of dr 
over all orientations is d2 d Vd.r d ~ 1 dr, Vd = (7r/4) d / 2 /r(l + 
d/2) being the volume of a (/-dimensional sphere of unit 
diameter. Note that in Eq. (|4.4[) the integration variable 

has changed from £ to a^2 and we have simplified the 
notation as aia — > coa and yla — > Uoa- This change 
avoids the need to specify the ^-dependence of <Jvl ■ 

Equation (|4.4[) applies to any choice of the set {a aj }. 
Henceforth we will assume the condition a ai > h(a a + 
<Jj), where a a and a 1 are the diameters of particles of 
species a and 7, respectively. This implies that the HS 
mixture is either additive or has a positive nonadditivity. 
In such a case, Eq. (|4.4j) can be further simplified. To 
that end, it is convenient to decompose /3^ x into two 
pieces, 



,CX.II 



(4.5) 



In the first contribution, the condition <ro a <\a a implies 
that the solute can lie "inside" a particle of species a. 
This makes the contribution Z?^™' 1 very easy to evaluate. 
First, by reversing the steps leading from Eq. (|2.10l) to 
Eq. (|4.4p . one can write 



(4-8) 



where £0 denotes the common value of the coupling pa- 
rameter £ at which <r£a = coa takes the value \o~ a for all 
a. According to Eqs. (|2~T4"]) . ([2~T5]l . and 



x f dr n e (l r «>- r j|-^iV 9 ) 

j=i 

Next, thanks to the condition er Q7 > ^(a a + cr 7 ), we note 
that in Eq. (|4.9|) the spatial regions excluded to particle 
i = by the remaining particles j — 1,...,N do not 
overlap, so that 

J dr [] e (ko - r, l - ^r ej ) = F(l - 77), (4.10) 

where n = Vd J2a Pa a a ^ s the total packing fraction of the 
solvent system. Thus, Eq. (|43| reduces to Qn+i(A ^0 = 
Q N (/3, ^)(1 - ri) and therefore = - ln(l - 77). 

Taking all of this into account, we finally get 



/3 f il x = - \n(l-Ti)+d2 d v d dcr 0Q <J^ a 1 y 0a (o- 0a )- 

(4-11) 



V. APPLICATION TO HARD-SPHERE 
MIXTURES IN THE PY AND SPT 
APPROXIMATIONS 



In this section, we apply Eq. (|4.11[) to the derivation of 
the chemical potential of a three-dimensional additive HS 
fluid mixture, i.e., a a -y = \{c a + er 7 ), as resulting from 
the PY and SPT approximations. Given an s-component 
mixture, the corresponding contact values are 



1 | 3 r] a a a 1 M 2 
1 - 77 2 (1 - r)Y a aj M 3 



+q 



o- a a 1 M 2 
(1 - nf 



(5.1) 
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where 



case unless q = |. This is because from Eq. (|5.5[) we get 



M n = 



(5.2) 



and the parameter q takes the values q = and q = | 
for the PY [13 and SPT @, & EMI theories, respec- 
tively. The more accurate Boubhk-Grundke-Henderson- 
Lee-Levesque (BGHLL) [l7l - [l9j expression corresponds 
to the intermediate value q = \. 



Now, in order to apply Eq. (|4.11[) . we assume that a 
solute particle i = is introduced in such a way that it 
interacts with a particle of species a through a HS poten- 
tial of range <xoa , the contact value of the corresponding 
cavity function being yoat^oa)- hi principle, we are free 
to choose <jQ a within the interval \o a < <Jo a < a va . On 
the other hand, if we want to make use of the approxima- 
tion (|5.1[) . we need to make the solute particle interact 
additively with the rest of the particles in the system. 
This is achieved by assuming that the solute particle is a 
sphere of diameter <jq within the range < (To < and 
coa = \ (on + (J a)- In that case, Eq. (|4. becomes 

PUT = - ln(l - rj) + — V X a / d<T CToa2/Oa((TOQ ! ), 

(5.3) 

where 7/oa(coa) is simply given by Eq. (|5.1j) with cr 7 — > 
ctq. After simple algebra, Eq. (|5.3|) yields 



tt 3 C: 3 



V 9 pJ Wq} U/W {p ^ } 24(1-77)3 — 



2 ^-2 2 



X (<7 a - er„) 1 



Aq 



(5.6) 



whereas Eq. (jTTTJ) implies the (Maxwell) symmetry rela- 
tion (|1.13[) . Therefore, except in the one-component case 
(a a = er„), dp v /dp a dn a ldp v unless q=\. Therefore, 
the PY prescription (|5.ip (with q = 0), when inserted 
into the chemical-potential route (|4.1ip , gives an expres- 
sion for the chemical potential of the mixture that is not 
strictly consistent with a well-defined free energy. On the 
other hand, this difficulty can be circumvented by using 
Eq. (HH) instead of Eq. $T7} to obtain [d{/3pa e *)/dp] {xa} . 
Integration over density (with the integration constant 
fixed by the condition lim p ^o flCX = 0) then yields 



f3a° 



ln(l - 77) + 



377 M 1 M 2 



3 ?7 2 



Mf 



1 - 77 M 3 



2(1 - r,f Mf 



6 - 977 + 2t7 2 | 6 ln(l-77) 
(1 — r]) 2 77 



1 4g\ 3M| 
3 J 2M| 



If now the chemical potential is obtained from Eq. (|5.7[) 
via Eq. (|1.7j) . the result differs from Eq. (|5.4j) unless g = 



(5.7) 



|, namely 



Pit 



= _ ln(1 _ ??) + ___|^ 



Mi 
M2 



37/ 



M2 
M 3 
Mfl 



77 Mi 



2(l- ?/ ) 
4g?7 2 



3M 2 1 - ?, M 3 3(1 - 7,) 2 



(5.4) 



Equation (|5.4[) shows that /3/i£ x is an explicit function 
of a v and depends on the diameters of all the species and 
on the partial densities through 77, Mi, M2, and M3, or 
cquivalcntly, through = ^ a /9 a <7™ = pM n with n = 
3. In terms of the latter quantities, Eq. (|5 -4[) becomes 



^ x = -ln(l-77) 



7T C2 



2 1 



T Cl 



2 1 



C: 



8 (1 - 77)2 



7T Co 



ClC2 



6 1 — 77 



+<1 



c 



54 (1 - 7 ? ) 3 



12 (1 - t?) 2 

(5.5) 



where one must take into account that 77 = ~£ 3 . The 
next step would be to derive the excess free energy per 
particle, a ox , from application of Eq. (|1.7p . Nevertheless, 
it turns out that the free energy associated with Eqs. 
(|5.4[) or (|5.5[) is not well defined in the multicomponent 



1 4g\ 9M| 
3 J 2M| 



6-9T7 + 277 2 6 ln(l - r?) 



M 2 3 
Ms' 1 ' 



(5.8) 



one 



Of course, although (dpa ox /dp^)p { Pa , u y 
has {dpa e */dp u ) 0{pa ^ } = E^/C^ so that Eq. 

(|5.7|) is recovered from both alternative expressions for 
the chemical potential via Eq. (|1.9p . Obviously, the right- 
hand side of Eq. (|5.8p satisfies the symmetry condition 
d 2 pa ex /dp a dp„ = d 2 pa ex /dp l/ dp a . 

The equation of state corresponding to the free energy 
(|5.7p is obtained from Eq. (|1.3p as 



1 



377 M X M 2 



37, 2 



M 3 



l-r? (I-77) 2 M 3 



(1 - 77)3 M| 



- 1- 



4g\ 3M 2 3 



2M 3 2 



6-I577+II7/ 2 ln(l-7i) 



(1 - 77)= 



6- 



(5.9) 



where we have taken into account that the ideal-gas con- 
tribution is linip^o Z = 1. In Eq. (|5.9p the subscript /i 
has been introduced to emphasize that this equation of 
state has been derived from the chemical-potential route. 
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Let us contrast Eq. (|5.9|) with the equation of state 
obtained from the virial route, Eq. (jl.ip . In the case of 
a (three-dimensional) HS fluid, Eq. becomes 



477 v - 

1 + w 3 ^ Xa 



a, 7 



Insertion of Eq. (|5.ip gives 

1 3?y MiM 2 



(5.10) 



Z v — 



3t7 2 



Af 2 3 



1 - n (1 



1 



1 



-r,f 
3 



(1 



7?)3 Mf 



(5.11) 



where the subscript v refers to the use of the virial route. 
Comparison between Eqs. (|5.9[) and (|5.11j) shows that the 
chemical-potential and virial routes coincide in the SPT 
(q = |) but not in the PY (q = 0) or BGHLL (q = ±) 
approaches. Interestingly enough, the SPT equation of 
state coincides with that derived from the PY solution 
[l2| via the compressibility route (|1.5[) . 

To summarize, the three PY equations of state as ob- 
tained from the chemical-potential, virial, and compress- 
ibility routes are 



vpY- 



= Z^q = 0), 



Z v (q = 0), 



9= 4 



(5.12) 



(5.13) 



(5.14) 



respectively. 

The celebrated Boubh'k-Mansoori-Carnahan- 
Starling-Leland (BMCSL) equation [U [H is obtained 
as an interpolation between the PY-v and the PY-c 
equations of state with respective weights ^ and |, i.e., 

^BMCSL = -ZpY-v + g^PY-c 



= Z v q = 



1 



(5.15) 



Having a third PY equation of state, Eq. (|5.12[) , it seems 
natural to construct an alternative interpolation formula 



^PY-/xc 



aZ 



+ {l-a)Z_ 



(5.16) 



In the one-component case 0] , values of the interpolation 
parameter a ~ 0.4 were seen to provide a better equation 
of state than the standard Carnahan-Starling equation 
[2l| . In particular, the values a = | and a = were 
explicitly considered. 

In order to assess the performance of Eqs. (|5.12[) 
(|5.16[) against computer simulations for binary mixtures, 
we have chosen the highest packing fraction (77 = 0.49) 
and the two largest size disparities (o^/fi = 0.6 and 




FIG. 1. (Color online). Plot of the compressibility factor Z 
as a function of the mole fraction xi for a binary mixture 
with a packing fraction rj = 0.49 and a size ratio 02/ ai = 0.6 
(top panel) or 0-2/0-1 = 0.3 (bottom panel). The symbols 
are computer simulation values [22], while the lines stand for 
(from top to bottom) Eqs. ((5TT4)) . ifFTTB) with a = 0.37, (pTTol) , 
(|5.12p . and (|5.13[) . respectively. 



CT2/cri = 0.3) considered in Ref. |22|. The simulation and 
theoretical results are displayed in Fig. [T] It is observed 
that, as expected, Zpy- V underestimates the simulation 
values, while Zpy_ c overestimates them. The chemical- 
potential route ZpY-fi lies below the simulation data, but 
it exhibits a better behavior than the virial route Zpy- v - 
The weighted average between Zpy- v and ^py-c made in 
the construction of the BMCSL equation of state (|5.15[) 
does a very good job. A slightly better agreement is 
obtained from the weighted average between Zpy-^ and 



Zpy-c [cf. Eq. ([5T6|) ] with a = | or a 



(not shown) 



but a value a = 0.37 provides especially accurate results. 



VI. CONCLUSIONS 

In this paper we have revisited the problem on the 
derivation of the chemical potential of a fluid system at 
equilibrium from the knowledge of the pair correlation 
functions. The result, Eqs. (|2.23p or (|2.24j) . applies to a 
mixture with any number of components and generic in- 
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teraction potentials <j> ai (r) between particles of species a 
and 7. The obtention of the chemical potential of species 
v relies on the formal introduction of an extra particle 
(the solute) coupled to the rest of the particles of the 
system (the solvent) via a series of interaction potentials 
<^ia(r)j which depend on a certain charging parameter 
£ in such a way that at £ = the solute ignores the 
presence of the solvent particle, while at £ = 1 the so- 
lute becomes indistinguishable from a solvent particle of 
species v [see Eq. (|2.13|) ]. The choice of the protocol 
leading from <0a (r) = at £ = to <0l (r) = 0„ Q (r) at 
£ = 1 remains arbitrary. 

In contrast to the other three conventional routes to 
thermodynamics [see Eqs. (|1.2[) . and 1)1.50 ]. which 

only need the pair correlation functions of the solvent 
system, the chemical-potential route requires the solute- 
solvent pair correlation functions for every value of the 
coupling parameter £. This implies that, even in the 
one-component case, this fourth route makes use of the 
pair correlation function of a binary mixture, albeit one 
of the species (the solute) is present with a vanishing 
concentration. This inherent multicomponent character 
of the chemical-potential route can in principle hamper 
its practical implementation. 

As is well known, thermodynamic quantities obtained 
from a common approximate pair correlation function via 
independent routes are not necessarily consistent. In par- 
ticular, the thermodynamic relations (|1.10|) - (|1.12|) can 
be violated when the left-hand sides are evaluated from 
Eq. ()2.23l) and the right-hand sides are evaluated from 
Eqs. ()1.2I) . ()l.ll) . and (|1.5|) . respectively. Moreover, the 
chemical-potential route introduces extra sources of pos- 
sible thermodynamic inconsistencies. First, in the case of 
a mixture, the symmetry condition (|1.13|) may not be ful- 
filled. This is a consequence of the fact that, in contrast 
to the energy, virial, and compressibility routes, where 
global quantities are obtained, the chemical-potential 
route provides a quantity for each separate component 
of the mixture. On the other hand, Eq. ()1.13l) is trivially 
satisfied by one-component systems. A second source 
of inconsistency is much subtler. As said before, the 
chemical-potential route implies to load a new particle 
into the system by means of a coupling parameter £ and 
it is not guaranteed that the final result will be indepen- 
dent of the protocol followed in the loading process, as it 
should be. From that point of view, given an approximate 
theory, this route is expected to provide a whole class of 
results, rather than a unique one, even for one-component 
systems. This protocol-related thermodynamic inconsis- 



tency is indeed observed in the PY solution of Baxter's 
sticky-hard-sphere model [23|], as will be extensively an- 
alyzed in a forthcoming paper 24 1. 

The general scheme has been particularized to the pro- 
totype HS fluid mixture. In this case, since the inter- 
action potential Q7 (r) depends on a single parameter 
<7 Q7 , it is possible to eliminate the coupling parameter £ 
in favor of the solute-solvent interaction range <7o a , thus 
avoiding the protocol problem alluded to above. The re- 
sult, given by Eq. (|4.4|) . is valid for any HS mixture, both 
additive and nonadditive. Further progress can be made 
if er Q7 > \(u a + <Ty), i.e., if the HS mixture is either ad- 
ditive or has a positive nonadditivity. In such a case, the 
evaluation of the contribution to the chemical potential 
when < o§ a <\a a amounts to the trivial computation 
of the volume accessible to a point particle in a sea of 
non-overlapping solvent spheres. 

As a practical implementation of the chemical- 
potential route, we have considered the SPT and PY 
approximations for three-dimensional additive HS mix- 
tures. The SPT result is consistent with the symmetry 
condition ()1.13p . while the PY result is not. This rep- 
resents a neat example showing that the internal consis- 
tency condition 130 is not guaranteed by an approxi- 
mate RDF. Moreover, the SPT chemical potential is con- 
sistent with the SPT virial equation of state (which coin- 
cides with the PY compressibility equation of state) , but 
the PY chemical potential yields an equation of state 
that differs from the PY virial equation. In fact, the for- 
mer turns out to be more accurate than the latter, as 
comparison with computer simulations reveals. Interest- 
ingly, an interpolation between the PY compressibility 
and chemical-potential routes provides slightly better re- 
sults than the conventional interpolation between the PY 
compressibility and virial routes giving rise to the BM- 
CSL equation of state. 
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